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ABSTRACT 

We  show  that  the  frequency  bias  caused  by  distributed 
cavity  phase  has  a  strong  dependence  on  microwave 
power.  We  also  show  that  frequency  biases  associated 
with  microwave  leakage  have  distinct  signatures  in  their 
dependence  on  microwave  power  and  the  physical 
location  of  the  leakage  interaction  with  the  atom. 

1.  INTRODUCTION 

The  subject  of  frequency  shifts  in  atomic  frequency 
standards  caused  by  distributed  cavity  phase  and 
microwave  leakage  in  the  microwave  spectrum  goes 
back  to  the  earliest  days  of  the  thermal  beam  standards 
[1,  2],  and  has  been  the  subject  of  continuing  theoretical 
and  experimental  work  over  the  last  fifty  years.  [3-8]. 
Laser-cooled  fountain  frequency  standards  pose 
problems  with  respect  to  both  distributed  cavity  phase 
and  microwave  leakage  different  from  those  associated 
with  thermal  beam  standards.  This  is  due  both  to  the 
very  different  microwave  structure  used  in  fountains  as 
well  as  the  low  center-of-mass  velocity  and  very  narrow 
velocity  distribution,  which  allows  operation  at 
significantly  elevated  microwave  power  in  fountain 
standards. 

We  discuss  the  dependence  on  microwave  power  of 
the  frequency  biases  induced  in  an  atomic  fountain  by 
the  distributed  cavity  phase  (DCP)  and  microwave 
leakage.  The  entire  discussion  takes  place  within  the 
limit  of  small  detuning  because  the  velocity  of  the 
atoms,  typically  a  maximum  of  ~3  m/s,  results  in  a 
maximum  Doppler  shift  of  only  -100  Hz.  We  first 
discussed  the  power  sensitivity  of  the  distributed  cavity 
phase  shift  (DCPS)  in  [9],  and  subsequent  work  both  by 
our  group  [10]  and  the  Gibble  group  at  Penn  state 
confirm  the  existence  of  a  power  dependent  DCPS  [11]. 

The  analysis  of  the  phase  of  the  microwave  field 
within  the,  typically  TE0n,  microwave  cavities  used  in 
cold-atom  fountain  frequency  standards  has  been 
extensively  analyzed.  DeMarchi  provided  the  seminal 
contributions  of  proving  a  correspondence  between  the 
phase  and  power-flows  within  the  microwave  cavity, 
identifying  the  preferred  cavity  configuration  and 
modeling  the  phase  of  the  microwave  field  within  the 
cavity  [12-15].  Three-dimensional  analyses  of  the 
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microwave  field  within  the  cavity  have  been  completed 
by  several  authors  [16-20].  In  many  of  these  studies 
large  phase  excursions  are  predicted.  At  least  one 
author  has  used  these  phase  excursions  to  predict 

frequency  shifts  in  the  df/S  >10~16  range  as  a  result  of 

distributed  cavity  phase  [17],  larger  than  the  current 
value  assigned  in  the  error  budget  of  many  operational 
fountain  type  primary  frequency  standards.  We  note 
that  the  theory  presented  here  allows  one  to  test  that 
prediction;  we  have  done  so  and  found  the  prediction  of 
[17]  incorrect. 

We  begin  by  briefly  reviewing  the  properties  of  the 
microwave  field  within  a  resonant  TE0n  cavity.  We 
then  solve  the  Schrodinger  equation  for  a  two-level 
atom  in  the  case  where  the  field  within  the  microwave 
cavity  has  both  real  and  imaginary  parts,  which  is  the 
case  required  to  understand  both  the  distributed  cavity 
phase  frequency  shift  and  the  microwave  leakage 
frequency  shift.  Finally  we  examine  the  dependence  on 
microwave  field  amplitude  of  the  frequency  bias  caused 
either  by  distributed  cavity  phase  or  microwave  leakage. 

2.  MICROWAVE  CAVITIES  AND  PHASE 

The  typical  microwave  cavity  used  in  cesium  (and 
rubidium)  fountain  frequency  standards  is  a  cylindrical 
cavity  resonating  in  the  TE0n  mode  at  the  hyperfine 
frequency  of  the  atom,  9.193  GHz  in  the  case  of  cesium. 
See  Fig.  1.  The  “z-axis”  of  the  cavity  is  aligned  with 
the  gradient  of  the  gravitational  potential  and  atoms 
enter  and  leave  the  cavity  via  below-cutoff  waveguides. 
This  cavity  has  the  crucial  property  of  allowing 
relatively  large  diameter  (2ra  -  1  cm  at  v0  =  9.2  GHz.) 
cylindrical  waveguides  for  atoms  to  enter  and  leave  the 
cavity  without  unduly  influencing  the  TE01 !  mode  of  the 
cavity  and  thereby  causing  large  phase  gradients  in  the 
microwave  field  [13].  To  lowest  order  the  field  within 
the  cavity  is  describable  as  purely  TE  ( Ez  =  0 )  and  all 
field  components  are  derivable  from  the  longitudinal 
magnetic  field  Hz(x,y)  [13,14].  Under  these 
assumptions,  and  using  a  trivial  extension  of  the 
notation  of  DeMarchi  in  [14],  the  longitudinal  field  can 
be  written  as 
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Figure  1.  A  schematic  diagram  of  the  TEou  microwave  cavity  used  in 
most  cesium  fountain  frequency  standards. 
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where  $?(x,y)  is  the  distributed  cavity  phase  under 

discussion  here.  The  expression  for  the  real  part  of  the 
field  within  the  cavity  is  given  (with  a  change  to 
cylindrical  coordinates  p.  <j)  and  z )  to  lowest  order  by 
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where  Jn  is  the  nth  Bessel  function  and  p' m  is  defined  as 


the  m  solution  of 


dJ„_(x) 

dx 


=  0 ,  rc  is  the  radius  of  the 


cavity  and  d  is  the  height  of  the  cavity  (see  Fig.  1  and 
[21]).  The  first  term  is  the  desired  TEon  fundamental 
mode  of  the  cavity,  while  the  infinite  sum  of  higher 
order  modes  with  coefficients  An>mj  comes  from  the 
required  cavity  imperfections,  especially  the  apertures 
on  the  end  caps,  which  allow  atoms  to  enter  and  exit  the 
cavity,  and  the  feed  apertures  or  antennas  that  couple 
power  into  the  cavity.  The  cavity  phase  is  defined  by 
the  relation 

^m(Hzp 
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Re(Hz) 
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where  the  real  part  of  H  is  approximated  by  Eq  (2),  or  to 
a  reasonable  approximation  that  allows  for  an  analytic 
result  the  real  part  of  H  as  approximated  by  the  first 
term  in  Eq.  (2).  The  imaginary  part  of  the  field  is  a 
result  of  losses  within  the  cavity,  primarily  ohmic  losses 
of  the  surface  currents  on  the  cavity  walls.  Whatever 
the  details  of  the  imaginary  part  of  the  field,  it  can  be 
expanded  as  an  infinite  sum  of  TEnml  modes  that 
provide  a  set  of  basis  functions  [19,21]. 
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where  8  is  the  skin  depth  inside  the  cavity.  Higher  order 
azimuthal  (« >  0)  modes  are  excited  by  the  feed 

structure  of  the  cavity,  while  higher  order  radial  (m  >  l) 

modes  and  higher  order  longitudinal  modes  (/  >l)  are 

excited  by  the  cavity  apertures  as  well  as  construction 
imperfections  [21].  The  coefficients  Anml  and  Bnml  for 
1=2,4,... 2n  should  be  quite  small  for  a  cavity  with  mid¬ 
plane  feeds;  however,  this  may  not  be  true  for  the 
“improved”  cavity  design  described  in  [20]. 
Additionally  with  mid-plane  fed  cavities  the 
predominant  power  flow  caused  by  coupling  imbalance 
is  orthogonal  to  the  atomic  trajectory,  resulting  in  a  high 
degree  of  rejection  of  the  first  order  Doppler  shift;  this 
feature  is  lost  in  cavities  not  using  the  midplane  feed 
such  as  those  in  [20],  which  exhibit  a  high  sensitivity  to 
first  order  Doppler  shifts. 

Several  things  can  be  seen  by  inspection  of  the 
previous  equations.  First,  the  phase  of  the  microwave 
field  within  the  cavity  is  independent  of  microwave 
power  (as  it  should  be).  Second,  the  microwave  phase 
can  become  relatively  large  when  the  real  part  of  Hz  is 
sufficiently  small,  or  equivalently,  the  Rabi  frequency 
of  the  atom  is  sufficiently  small.  What  is  ultimately  of 
interest,  of  course,  is  not  the  value  of  the  phase  angle  of 
the  microwave  field,  but  rather  the  value  of  the  phase 
angle  imposed  on  the  atomic  wave  function  as  a  result. 
The  effect  on  the  atom  of  the  imaginary  part  of  Hz  is 
dependent  on  the  microwave  power,  in  spite  of  the 
phase  of  the  microwave  field  being  independent  of  that 
same  power. 

We  have  previously  expanded  the  imaginary 
part  of  the  field  using  a  power  series  [22,15],  which 
results  in  an  expression  such  as 


where  again  the  odd  terms  in  the  series  are  small  as  a 
result  of  symmetry  considerations.  The  Penn  State 
group  has  obtained  similar  but  more  complete  power 
series  expansions  [20].  The  expansion  of  the  field  in 
terms  of  a  normal  mode  picture  (Eq.  (4))  or  a  power 
series  (Eq.  (5))  yields  similar  power  sensitivities  of  the 
frequency  bias. 

In  order  to  quantify  the  effect  of  the  imaginary 
part  of  the  microwave  field  upon  the  atom,  we  now 
obtain  a  solution  to  the  time -dependent  Schrodinger 
equation  as  the  atoms  pass  through  a  cavity  with  fields 
described  by  Eq.  (2),  (4)  &  (5). 
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3.  SCHRODINGER  EQUATION  AND  RAMSEY 
LINESHAPES 

We  extensively  employ  the  theoretical 
framework  developed  by  Shirley  et  al.  [23-25]  and 
present  here  the  extensions  required  to  handle  both  the 
real  and  imaginary  phases  of  microwave  field. 

The  Hamiltonian  for  the  system  can  be  written  as  (cf. 
(7)  of  [25]) 


f  CO i  2bcos<yt  +  2b'sinfflt) 

n=h\  ,  ,  (6) 

^2bcosftlt  +  2b  sin®t  COb  J 

where  fr(Oa  and  ftcob  are  respectively  the  energies  of  the 
upper  and  lower  states.  The  interaction  Rabi  frequency 
for  the  real  part  of  the  microwave  field  is  given  by 
2b  =//BgRe(Hz)//i,  where  |iB  is  the  Bohr  magneton,  g 

the  Lande  g-factor  and  Hz  the  microwave  magnetic  field 
parallel  to  the  quantization  axis  imposed  by  the  external 
c-field.  A  similar  expression  applies  for 

b\  2b'  =  //BgIm(Hz)/?i,  the  Rabi  frequency  due  to  the 
imaginary  part  of  the  microwave  field.  Both 

b(t)  andb'(t)are  time -dependent,  owing  to  the  atoms 
motion  in  the  cavity. 

In  the  rotating  wave  approximation,  the 
Hamiltonian  in  Eq.  (6)  is  written 

f  coa  (b  +  ibV^l 

*  =  d  }  <7> 

Note  the  sign  change  in  b'(t)  in  the  off-diagonal 
couplings.  This  comes  about  because  the  rotating  wave 
approximation  selects  one  exponential  from  sin  cot  in 
one  coupling  and  the  other  exponential  in  the  other 
coupling  (Compare  to  (7)  and  (8)  in  [21]).  Using  the 
“phase  factored”  solutions,  a  and  [3,  (cf  (9  &  10)  of 
[24])  gives  us,  finally,  the  Schrodinger  equation  for  the 
system, 
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With  the  initial  condition  or(0)  =  1  and  /?(0)  =  0,  a  is 
the  probability  amplitude  that  the  system  remains  in  its 
initial  state  and  [3  the  probability  amplitude  that  the 
system  changes  state.  A  is  half  the  detuning  8(0  from 

,  1  ,  ,8co 


atomic  resonance 


where  ty0  =  (0^  -  (Oh  is  the  hyperfine  splitting  of  the 


atom.  A,  b,  and  b'are  all  real,  possibly  time-dependent, 
quantities. 

We  have  given  a  solution  to  (8),  valid  through 
first  order  in  A  and  b',  under  the  assumption  that  the 
detuning  is  small  compared  to  the  Rabi  frequency,  in 
[9].  We  quote  here  only  the  final  results. 

Using  our  solutions  mentioned  above  and  assuming  that 
the  real  part  of  the  excitation  is  the  same  for  both 
Ramsey  pulses  leads  to  a  transition  probability  of 


bo  is  defined  by  b0z  =  ^b[t)dt ,  that  is,  b0  is  the  average 

0 

value  of  Z>(t)  during  the  excitation  time  T,  which 


amounts  to  at  optimum  power,  and  g,  2  and  rji  2  are 

proportional  to  the  imaginary  part  of  the  microwave 
field  b';  the  subscripts  refer  to  the  first  and  second 
passage  through  the  Ramsey  cavity.  That  is, 


From  Eq.  (9)  the  transition  probability,  P,  is 
clearly  a  normal  Ramsey  fringe  (the  cos  ScoTR  term) 
plus  an  underlying  fringe  of  small  amplitude  and 
n!  2  displacement  (the  sin  8(oTR  term).  The  latter 
asymmetrically  distorts  the  Ramsey  curve  and  leads  to  a 
frequency  bias.  The  frequency  bias  is  proportional  to 
the  difference  between  the  transition  probabilities  PL 
and  PR  on  the  left  and  right  sides  of  the  central  Ramsey 
fringe  at  equal  detunings.  This  difference  is  given  by 


PL-PR  =  sin2  (2b0r) 


g2-g,  t  i)2  +  nt 
sin(b0r)  cos(b0r) 


(11) 


where  PL  R  is  the  probability  given  by  (9)  on  the  left  and 
right  sides  of  the  central  Ramsey  fringe,  respectively 

(that  is  with  8co  =  +  ^T  respectively).  Eq  (11)  is 

plotted  below  in  Fig  2.  The  complicated  nature  of  the 
power  dependence  of  the  frequency  bias  associated  with 
distributed  cavity  phase  is  immediately  apparent  in  Fig. 
2.  This  signature  should  allow  measurement  of  the 
frequency  bias  associated  with  distributed  cavity  phase, 
or  at  least  allow  the  placing  of  an  upper  limit  on  the 
effect. 

Two  points  should  be  stressed  here:  First,  as 
noted  in  our  earlier  work  [9],  a  power  dependence  of  the 
DCPS  induced  frequency  bias  occurs  whenever  the 
curvature  of  the  real  and  imaginary  parts  of  the  field  are 
different  over  the  atomic  trajectory.  This  is  bound  to 
occur  in  any  real  fountain  cavity,  which  must  have 
apertures  in  the  end  caps  for  the  atoms  to  enter  and 
leave.  Second,  a  frequency  shift  can  only  be  introduced 
by  a  field  in  quadrature  to  the  principal  part  of  the  field 
(that  is,  the  imaginary  part  of  the  field).  Variations  in 
the  real  part  of  the  field  cannot  lead  to  a  frequency  shift 
in  and  of  themselves. 
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Microwave  Excitation  (units  are  b0x=7t/2) 

Figure  2.  The  distortion  of  the  Ramsey  fringe  plotted  caused  by 
distributed  cavity  phase  with  respect  to  the  microwave  excitation 
(number  of  /r/2  pulses).  The  distortion  is  proportional  to  the 
frequency  bias  and  the  data  shown  here  are  normalized  to  the  shift  at 
optimum  power  (one  71  j 2  pulse). 


4.  MICROWAVE  LEAKAGE 

Microwave  leakage  induced  frequency  biases  have 
previously  been  investigated  by  several  authors  [6,7]. 
Because  the  earlier  work  was  primarily  in  the  context  of 
thermal  beam  clocks  operation  at  optimum  power  was 
investigated  and  the  power  sensitivity  of  the  effect  was 
not  elucidated.  Here,  because  we  are  primarily 
interested  in  the  effects  of  microwave  leakage  in  the 
context  of  fountain  frequency  standards,  we  investigate 
the  effect  in  the  context  of  small  detuning  (small 
Doppler  shift)  and  large  microwave  power. 

We  can  use  the  Hamiltonian,  Eq.  (7),  also  to 
investigate  microwave  leakage,  except  that  in  this  case 
both  b  and  b'  are  presumably  small.  We  divide  our 
analysis  into  three  cases:  leakage  before  the  two 
Ramsey  interactions,  leakage  between  the  Ramsey 
pulses,  and  leakage  after  the  Ramsey  interactions. 
Because  the  state-selection  mechanism  in  NIST-F1  and 
IEN-CSF1  destroys  any  microwave  coherence  and 
projects  the  atoms  into  a  pure  F=3  state  we  can  ignore 
the  case  of  leakage  before  the  first  Ramsey  interaction. 
We  note  however  that  the  solution  to  the  case  of  leakage 
before  the  first  Ramsey  interaction  can  be  obtained  from 
the  case  of  leakage  after  the  second  Ramsey  interaction 
by  the  substitution  of  —t  for  t. 

Some  general  results  can  be  obtained  [26]: 

1.  Microwave  leakage  in  phase  (the  b  term)  with 
the  field  in  the  Ramsey  cavity  does  not  (indeed 
cannot)  cause  a  frequency  shift. 

2.  Microwave  leakage  in  quadrature  phase  (the 
b'  term)  does  not  cause  a  frequency  shift  if  it  is 
applied  symmetrically  with  respect  to  the 
center  of  the  Ramsey  interaction  at  TR/2. 


3.  Microwave  leakage  in  quadrature  phase 
applied  asymmetrically  with  respect  to  TR/2 
excites  a  Ramsey  fringe  shifted  with  respect  to 

the  central  fringe  (asin^ryT), )  term)  much  like 

that  in  Eq.  (9),  which  can  cause  a  frequency 
bias. 

We  investigate  the  power  dependence  of  these 
frequency  biases  next.  Full  solutions  to  the  Schrodinger 
Equation  for  all  three  leakage  cases  are  given  in  [26]. 
Below  the  various  specific  cases  of  microwave  leakage 
are  examined. 
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Ramsey  Excitation  -  in  n/2  pulses 


Figure  3.  The  modeled  frequency  shift  induced  by  leakage  above  the 
Ramsey  cavity  in  NIST-F1  and  IEN-CSF1.  The  solid  curve  is  the  case 
where  the  leakage  amplitude  is  constant  while  the  field  in  the  Ramsey 
cavity  is  increased.  The  broken  curve  is  the  case  where  the  leakage 
field  amplitude  is  proportional  to  the  amplitude  of  the  field  within  the 
Ramsey  cavity.  The  microwave  leakage  power  of  each  curve  was 
chosen  for  convenience  and  has  no  particular  significance. 


5.  MICROWAVE  LEAKAGE  BETWEEN  THE 
TWO  RAMSEY  PULSES 


The  dependence  of  the  frequency  bias  on  the 
amplitude  of  the  microwave  field  in  the  Ramsey  cavity 
is  shown  in  Fig  3.  It  is  immediately  apparent  that  the 
frequency  bias  caused  by  microwave  leakage  is  zero  at 
integer  multiples  of  the  optimum  excitation  amplitude, 

b0T=  (2^  +  1)^/  ^  n=0,l,2...  .  The  residual  frequency 

bias  can  be  measured  with  some  leverage  by  measuring 
well  away  from  optimum  power,  as  illustrated  in  Fig.  4. 
In  NIST-F1  the  microwave  power  is  set  to  within  ±0.1 
dB  of  optimum,  denoted  by  the  dotted  lines  in  Fig.  4. 
The  frequency  shift  induced  by  microwave  leakage  at  ± 
3  dB  from  optimum  excitation  is  some  30  times  greater 
than  the  shift  expected  within  0.1  dB  of  optimum.  The 
measurement  at  ±  3  dB  therefore  gives  a  “leverage”  of 
30  over  the  standard  operating  conditions,  at  least  for 
the  specific  case  of  leakage  above  the  Ramsey  cavity. 
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Figure  4.  The  frequency  bias  from  microwave  leakage  above  the 
Ramsey  cavity  can  be  measured  in  a  “leveraged”  fashion  by 
operating  well  away  from  optimum  microwave  power  in  the  Ramsey 
cavity.  Assuming  that  the  normal  operating  conditions  keep  the 
microwave  power  within  0.1  dB  of  optimum,  a  leverage  of  about  30 
can  be  obtained  by  operating  3  dB  above  or  below  optimum  instead. 

6.  MICROWAVE  LEAKAGE  AFTER  THE 
RAMSEY  INTERACTION 

In  the  case  of  both  NIST-F1  and  IEN-CSF1 
this  is  the  most  likely  location  for  microwave  leakage  to 
interact  with  cesium  atoms.  The  power  dependence  of 
this  shift  is  shown  in  Fig.  5.  The  behavior  of  this  shift  is 
distinct  from  that  caused  by  leakage  above  the  Ramsey 
cavity  in  that  its  frequency  with  respect  to  the  Rabi 
frequency  is  half  that  of  the  case  of  leakage  above  the 
Ramsey  cavity.  This  signature  can  be  used  to  identify 
the  source  of  a  frequency  bias  as  being  caused  by 
leakage  either  above  or  below  the  Ramsey  cavity. 
Unfortunately  this  shift  is  difficult  to  distinguish  from 
the  frequency  shift  caused  by  distributed  cavity  phase, 
as  shown  in  Fig  6.  Various  approaches  to  separate  these 
two  effects  are  discussed  in  [26]. 


Ramsey  Pulses  (units  of  b0x=7t/2) 


Figure  5.  The  modeled  frequency  bias  induced  by  leakage  below  the 
Ramsey  cavity.  The  solid  curve  is  the  case  where  the  leakage 
amplitude  is  constant  while  the  field  in  the  Ramsey  cavity  is  increased. 
The  broken  curve  is  the  case  where  the  leakage  field  amplitude  is 
proportional  to  the  amplitude  of  the  field  within  the  Ramsey  cavity.  In 
both  cases  the  microwave  power  level  has  been  chosen  purely  for 
convenience  and  has  no  particular  significance. 


The  frequency  biases  caused  by  distributed  cavity 
phase  and  microwave  leakage  have  distinct  signatures 
when  measured  as  a  function  of  the  amplitude  of  the 
microwave  excitation  in  the  Ramsey  cavity.  The 
various  shifts  can  be  measured  and  either  “fixed”  by 
eliminating  the  source  of  the  leakage  (best),  or  perhaps 
calibrated  and  corrected. 
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Figure  6.  A  comparison  of  the  power  dependence  of  the  shift  caused 
by  microwave  leakage  (dotted  curve)  and  the  shift  caused  by 
microwave  leakage  below  the  Ramsey  cavity  (solid  curve)  as  a 
function  of  the  microwave  excitation  of  the  atoms  in  the  Ramsey 
cavity. 

It  has  been  assumed  here  that  the  leakage  terms  can  be 
treated  in  the  small  detuning  limit.  Various  other 
approaches  such  as  using  the  nonresonant  AC  Stark 
shift  have  been  applied  to  the  case  of  microwave 
leakage  in  a  cesium  fountain  [27].  We  hold  that  these 
cases  are  not  valid  approximations  in  a  cesium  fountain; 
the  Doppler  shift  caused  by  atomic  motion  is,  at  100  Hz 
maximum,  much  too  small  to  be  considered  within  the 
framework  of  an  off-resonant  AC  Stark  shift.  The 
frequency  shift  caused  by  microwave  leakage  in  a 
fountain  depends  on  the  amplitude  of  the  microwave 
field  and  not  on  the  power  (square  of  the  amplitude). 
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